Anomalous resonance phenomena of solitary waves with internal modes 
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We investigate the non-parametric, pure ac driven dynamics of nonlinear Klein-Gordon solitary 
waves having an internal mode of frequency fli. We show that the strongest resonance arises when 
the driving frequency 5 — Q.i/2, whereas when 5 = Q,i the resonance is weaker, disappearing for 
nonzero damping. At resonance, the dynamics of the kink center of mass becomes chaotic. As we 
identify the resonance mechanism as an indirect coupling to the internal mode due to its symmetry, 
we expect similar results for other systems. 
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An important paradigm established over the last two 
decades is that solitary waves or solitons behave very 
much like point particles when subjected to (a large class 
of) external forces and perturbations Q-D- However, 
many solitary waves possess one (sometimes, more than 
one) internal or shape mode |^J^, and in that case the 
particle picture of their dynamics may be oversimplified: 
Indeed, the internal mode can temporarily store energy 
and release it at a later stage, giving rise to resonance 
phenomena in solitary wave collisions pi or in solitary 
wave interactions with inhomogeneities ^ . As internal 
modes are quite common in nonlinear systems, either in- 
trinsically or a result of small perturbations [Q , the ques- 
tion of their influence on the dynamics of solitary waves 
is a very general and relevant one. 

One aspect of solitary wave dynamics that has proven 
itself difficult to understand is that of topological solitary 
waves or kinks subjected to pure, i.e., non-parametric 
ac driving. Only recently |^ the ac driven dynamics of 
sine-Gordon kinks has been definitely clarified, settling 
a question that had received different and contradictory 
answers. However, since the sine-Gordon kink does not 
have an internal mode, it does not provide information 
about the interplay of those modes with external ac driv- 
ings, which, to our knowledge, has only been analyzed 
for breathers, non-topological excitations with an inter- 
nal mode (see |^ and references therein). Naively, the 
only new phenomenon one expects from a non-parametric 
external driving is a resonance when its frequency, 5, 
matches that of the internal mode, fl^. The aim of this 
letter is to show that, in fact, the actual scenario is most 
unexpected and highly non-trivial. As we will see below, 
a strong, anomalous resonance arises when 5 — f2i/2, 
whereas the normal resonance at S — fli is definitely 
weaker, only possible at exactly zero damping, and even 
then it can be suppressed by appropriate choices of other 
parameters. We expect this result to be generic, because 



our analytical approach allows us to identify the mecha- 
nism for such a peculiar phenomenon: The ac force does 
not act directly on the internal mode (because of symme- 
try reasons), but rather, they interact indirectly via the 
translational motion which couples to the internal mode. 
Our predictions are fully confirmed by numerical simu- 
lations, which in addition show the implications of these 
resonances for the kink dynamics. 

As a specific example of a kink with internal mode, we 
take the well known 0^ equation, which, when driven 
with an ac force f{t) = €sm{6t + 6q), reads 

c^tt-<i^.. + U'{c^)^-(icj)t + I{t), (1) 

where U{4)) — (0^ — 1)^/4 and /3 is a damping coefhcient. 
Previous related works on this system are ||^ , where reso- 
nances in the presence of an external (time independent) 
potential have been considered, and jl^, that dealt with 
non-resonant, high frequency parametric ac drivings. For 
our problem, our analytical approach will be the well 
known collective coordinate (CC) method [Hj^. A first 
order of approximation is given by the McLaughlin-Scott 
method : We assume that the solution of (|l|) is of the 
form 

i:-X{t) 



{x, t) = tanh 



(2) 



-io^Ji-v{tY 

where Iq = \/2- The center of the kink X{t) and its ve- 
locity V{t) are related by X{t) = /J dt'V(t') + X{0), and 
both are unknown functions describing the motion of the 
kink as a coherent entity. By means of a standard proce- 
dure involving conservation laws, an ordinary 
differential equation of motion for V{t) can be obtained, 
linearized and solved, yielding 

r{t) 



V{t) 



r{i) 
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X [/? sin{St + Sq)-S cos(5t + Sq)] 
3V2e 



c = 7ol/(0) + 



Psm{6o) - 6cos{6o)l (4) 



where 70 = l/i/l — ^(O)^. For the undamped case 
{(3 — 0), we see that the kmk oscihates if 'yoV{Q) — 
3\/2e cos((5o)/(2(5); otherwise, the kink will move either 
to the right or to the left, depending on the relation be- 
tween the parameters of the ac force and the initial ve- 
locity. Note also that such a dc motion is absent when 
(3^0. These results were numerically confirmed for sine- 
Gordon kinks (which do not have internal modes) in . 

In order to include internal mode effects, we proceed 
as follows: We rewrite Eq. (m as 



^ = -^-P^ + fit), 



5H_ 



where ip = cj), the dot meaning derivative with respect to 
time, and 



H = 



f 1 1 -1 



(6) 



is the Hamiltonian of the system when e — [3 — 0. We 
now make the Ansatz (j){x, t) — (f)[x — X{t), l{t)], whereas 
from the definition of ^ we have that 'ip{x,t) = ip[x — 
X{t),l{t),Xj]. As in the McLaughlin-Scott method, 
X{t) represents the kink center position, but now we in- 
troduce a second collective variable l{t), that will stand 
for the kink width, i.e., the internal mode excitation, be- 
low. 

The procedure to obtain the CC equations correspond- 
ing to this generalized travelling wave Ansatz has been 
put forward in Basically, it consists of insert- 

ing our Ansatz into (|5|), multiplying the first equation 
by d(j)/dX and the second one by d4>/dX {d(j)/dl and 
dijj/dl), taking their difference and integrating over x. 
This yields a rather cumbersome, ordinary differential 
equation for X(t) [l{t)], which we omit here for brevity. 
The next step is to choose a specific functional form for 
(j), which we do following the work of Rice p3| , and let 



{x~X(t),l{t)]^<j>o 



X{t) 



lit) 



(7) 



where (l)o[{x — X)/lo] is the static kink solution, which is 
and odd function (with respect to its center) for the (j>'^ 
and for any other even potential. Upon particularization 
of the CC equations for this form for (p we finally obtain 

X xi X 

Molo- - M^l^— = -pMoh- - 2f{t), (8) 
aMolo\ + Mok^ = K"'\l, U X) - paM^lJj- (9) 
where K"'* = —dE/dl and 



(10) 



is the kink energy. Importantly, in obtaining Eq. (^) a 
term of the form f(t) J_^dxd(j>Q/dl, coming from the 
coupling of the ac driving to the internal mode, has 
vanished because of symmetry. For the 0^ equation, 
(j)o = tanhx in Eq. (0), which yields a = (tt^ - 6)/12 
and Mo = 2^/2/3. 

Let us now simplify these expressions in order to make 
its physical significance more transparent. To begin with, 
Eq. (^) can be solved for X/l, yielding: 



p _ MqIqX _ ^^l3sm{5t + So) - dcos{St + Sq) 



2e 



(/32 + 52) 
^sin(<5o) - (5cos((5o) MoloV{0) 



{f3^ + <52) 



(11) 



where Is = Iq /70 • Inserting Eq. (O) into Eq. (M , we find 



a[f - 211 - 2f]ll] = -2 
'0 



Ml 



- 1. 



(12) 



Although Eqs. ( pT| ) and (|T^) are quite complicated, 
choosing 1^(0) and the phase 60 so that the exponential 
terms in Eq. (11) vanish, we can use a change of variables 
(proposed in |15] for the dc driving case) to transform the 
Eq. (|l^) into a Pinney-like equation (see [|l^ and refer- 
ences therein), which can be solved in terms of Mathieu 
functions ii (3 — (when (3 ^ we have not been 
able to solve this problem analytically). 

In any event, we do not need the analytical expressions 
for the solution of the system (11), ( p^ ) to understand 



the physics predicted by our approach. The term P in 
Eq. (12) is an oscillatory function with frequency 26 [see 
Eq. (11)]; hence, we can immediately expect a resonance 
when the external frequency 6 is half the frequency of the 
internal mode, ^ 1/^/aZo in the Rice approximation. 
For the cf)^ model VIr overestimates Vli = \/3/2 by 1.7% 
[0. The analytical solution for (3 = confirms this ex- 
pectation, while numerical integration of Eq. ( |l2|) proves 
that the behavior of I is that of a resonant, damped oscil- 
lator when (3 0. When 5 = Hr, inspection of Eqs. ( |Tl| ) 
and (12) leads to the conclusion that another resonance 
should be found at 6 = Qr only if (3 = Q (otherwise it is 
a transient phenomenon of lifetime (3~^)] even then, by 
choosing V^(0) and (5o to cancel the non-oscillatory terms 
in Eq. (^) the resonance is completely suppressed. Fur- 
thermore, both analytically and numerically we have ver- 
ified that, far away from the resonances, the behavior of 
the kink center, X{t), is practically the same as the one 
predicted by the McLaughlin-Scott approach, Eq. (|^). 
This means that, within the CC framework, the behav- 
ior of ac driven (f)^ kinks is described by the McLaughlin- 
Scott Ansatz, and only for drivings close to VLr/2 (and 
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fill if /? = 0) such approach fails and resonant phenomena 
arise. 

This far, the CC method has provided specific predic- 
tions about the system behavior, but in order to assess 
their physical reality, we must verify that they arise in 
the partial differential equation (PDE) , Eq. (|l]) . The rea- 
son for this is twofold: First, underlying CC methods is 
the assumption that no (or a negligible amount of) ra- 
diation is generated by the perturbation, an assumption 
whose validity can only be assessed through comparison 
with the corresponding PDE. Second, even if that is the 
case, the CC equations predict an unbounded growth of 
l(t) at resonance, and it is difficult to understand what 
does that mean in physical terms for the kink of the full 
PDE, whose width is controlled by the properties of the 
equation. In view of this, we have computed the numer- 
ical solution of the PDE d^) by using the conservative 
Strauss- Vazquez scheme with Aa; = 0.1, At = 0.01, 
and a total system length of L = 400. Our initial con- 
dition was a kink at rest and Sq = 7r/2, a pair of values 
for which we should not see a resonance at fli . We come 
back to this below. We monitored the position and the 
velocity of the kink center as well as the total energy in 
the system, computed from the Hamiltonian (^). We also 
tried to measure directly the kink width, but we found 
that it is quite complicated to estimate it from the nu- 
merics, this being the reason why we have resorted to less 
direct measurements. 



both magnitudes is periodic, whereas at the resonance it 
becomes chaotic. Specifically, the energy increases with 
time: The closer to the resonance value, the faster the 
increment. The center motion is initially periodic, until 
the internal mode amplitude has increased too much and 
stored too much energy, subsequently releasing it through 
its coupling with the translation mode (which we know 
that exists from ) , eventually yielding the kink motion 
erratic as this process is repeated once and again. This is 
a clear evidence in favor of the resonance predicted from 
the CC treatment. We have verified that at resonance 
the kink motion is very sensitive to changes in the initial 
conditions, hence our claim of the appearance of chaos. 
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FIG. 1. Results for the PDE (|lD with /3 = 0, e = 0.01. (a) 
Total energy when 5 — 0.6100 (upper line), 5 = 0.6080 (lower 
line), (b) Same for the kink center, X{t). Insets show the 
same for S — 0.6102, a value for which the kink reaches the 
boundary of the numerical system before t — 25 000. 

Figure ^ shows examples of the kink center dynamics 
and its energy evolution both close to and away from 
the predicted resonance. Off-resonance, the behavior of 



FIG. 2. Resuhs for the PDE (|l|). Mean value of the energy 
(see text for the way it is computed) vs driving frequency S, 
(a) close to Q.i/2, (b) close to In both cases, the upper 
line corresponds to /3 = 0, and the lower line to 13 = 0.001. 

Figure ^ depicts the resonance as seen through the 
mean energy, computed as a time average from t = 10 000 
(after transients have died out in the damped case) to the 
end of the run at t = 25 000. Figure ||(a) shows the be- 
havior of this magnitude around 51^/2. It is clear from 
the plot that there is a strong resonance dX 5 k, 0.6102 w 
rii/2, both with and without dissipation. In fact, some 
points are missing in the (3 = Q line because the corre- 
sponding kinks, in their apparently random motion, left 
the system before the end of the run. On the other hand, 
in spite of our choice of initial conditions, for which no 
resonance is predicted at for /? = 0, a weak resonance 
can be seen at (5 w 1.225 ~ Vli in Fig. ^b). We note 
that the other two small peaks are spurious, since they 
appear, disappear, or change location depending on the 
choice of the length for the numerical simulation. We be- 
lieve that this discrepancy might come from the difficulty 
of numerically tuning the condition for its suppression. 
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In any event, that would be a very special case, and the 
behavior we find for these parameters is representative of 
what occurs for other choices of the initial velocity and 
the phase. Those give rise to a similar behavior, while, 
remarkably, the resonance at ilj is always weaker and nar- 
rower than that at 0^/2 (cf. Fig. |^). We also see that the 
prediction that the resonance at f2i is suppressed by the 
dissipation is also confirmed, the small peak in the lower 
line stemming from the transients, which are not exactly 
zero for 10 000 < i < 25 000. Another interesting remark 
is that Fourier analysis shows that the monotonous in- 
creasing of the energy that appears in Fig. ||(b) comes 
from the fact that, when 5 ^ 1.1, the lowest phonon 
mode (with frequency ujp = -\/2) begins to be excited, 
the amplitude of its excitation monotonically increasing 
as S ^ V^. No evidence for lowest phonon mode exci- 
tation is seen for the resonance at 0^/2; therefore, this 
is indeed a phenomenon arising from the coupling of the 
translation and the internal mode as predicted by the CC 
calculation. 

In summary, we have studied how ac forces affect soli- 
tary waves of kink type with an internal mode of fre- 
quency Qi. Specifically, we have clearly shown that the 
behavior of (j)'^ kinks under ac driving is very well de- 
scribed by the two- variable CC theory we have developed 
here. The main feature of the non-parametrically, ac 
driven 0* kink dynamics is that the strongest resonance 
occurs at (5 = rii/2, and not at the frequency one would 
expect, S — Qi. We emphasize that this novel resonance 
phenomenon is totally unexpected from the knowledge of 
the internal mode frequency, and arises from the indirect 
interaction of the external force with the internal mode 
via the translational motion. Although our results have 
been obtained for a specific example, the (f>'^ equation 
(|I|), other models with even on-site potentials and inter- 
nal modes, such as the double sine-Gordon equation, for 
instance, will behave similarly because a CC approach 
will lead to analogous results. The resonance found at 
the CC level manifests itself at the PDE level as erratic 
or chaotic (strongly dependent on the initial conditions) 
motion of the kink as the kinetic energy of the center 
of mass is stored into, and recovered from, the internal 
mode. 

On closing, we note that the anomalous resonances de- 
scribed here are important on their own, as examples of 
the highly non-trivial behavior of nonlinear systems and 
as hints about the mechanisms governing kink dynamics. 
In addition, we think that this phenomenon should be 
very general, in view of the recent finding that per- 
turbations of kink-bearing nonlinear systems often lead 
to the development of an internal mode. As this occurs in 
the discrete sine-Gordon model [Q, a resonance like the 
one discussed here could be relevant for the mode locking 
phenomena reported for that system in ]l9[ |. Finally, by 
using this discreteness induced internal mode, we point 
out that the resonance we find could be observed in ex- 



periments by using Josephson junction arrays as in pO| . 
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